A numerical study of primordial magnetic field amplification by infiation-produced 

gravitational waves 
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We numerically study the interaction of inflation-produced magnetic fields with gravitational 
waves, both of which originate from quantum fluctuations during inflation. The resonance between 
the magnetic field perturbations and the gravitational waves has been suggested as a possible mech- 
anism for magnetic field amplification. However, some analytical studies suggest that the effect 
of the infiationary gravitational waves is too small to provide significant amplification. Our nu- 
merical study shows more clearly how the interaction affects the magnetic fields and confirms the 
weakness of the influence of the gravitational waves. We present an investigation based on the 
magnetohydrodynamic approximation and take into account the differences of the Alfven speed. 
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I. INTRODUCTION 



Recently, various observations have confirmed the pres- 
ence of magnetic fields both in individual galaxies and on 
even larger scales such as in clusters of galaxies [1] . These 
magnetic fields have a typical strength of 10~^~~^G. 
Such large and strong magnetic fields are expected to 
have been amplified to their present levels by galactic 
dynamo theory [2[ . However this efficiency is a matter of 
some debate, because this mechanism requires seed mag- 
netic fields greater than ~ 10~^°G 0] and existing theo- 
ries for generating seed fields do not achieve the required 
strength Q. To fill the gap betvifcen observations and 
the dynamo theory, a large number of theoretical mod- 
els have been suggested to generate stronger seed fields 
or to amplify the weak seed fields generated from exist- 
ing theory. The minimum amplitude necessary in order 
for dynamo theory to produce the fields observed at the 
present epoch is estimated as approximately 10~'^*G at 
a comoving length of lOkpc 

In this paper, we focus on the amplification of inflation- 
produced magnetic fields by infiation-produced gravita- 
tional waves; a mechanism that was originally suggested 
by Tsagas (2003) [6^. Providing no exotic physics is as- 
sumed, infiation generally predicts quite small magnetic 
fields of roughly 10^^'^G [7|. Various mechanisms to 
break the conformal invariance of electrodynamics have 
been proposed to generate larger magnetic fields, but 
these are not widely accepted since they require the in- 
clusion of new physics beyond the standard model. Grav- 
itational generation of the massive Z-boson field has been 
suggested as a mechanism to naturally break conformal 
invariance within the standard model Q. However, 
the strength of the generated magnetic field with a rea- 
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sonable model of field evolution is estimated to be only 
10~^^G on scales of lOOpc, which is not large enough to 
explain observations of the present day universe. Orig- 
inally, Ref. 6] reported that such small seed magnetic 
fields may be amplified to the order of 10^^ by amplifi- 
cation induced by gravitational waves. The key idea of 
this amplification mechanism is the resonance between 
gravitational waves and magnetic field waves as a result 
of the interaction between them, which occurs only when 
these have the same wavenumber. Since inflation gener- 
ates both gravitational waves and magnetic flelds at all 
scales, this mechanism is considered to be able to work 
on the inflationary magnetic fields. However, this work 
has since been revisited by several papers |9l-(l]| which 
conclude that it does not provide such significant ampli- 
fication. 

In contrast to the majority of these studies, which have 
been performed analytically, we here tackle the problem 
with a numerical approach. The advantage of such a nu- 
merical study is that it enables us to accurately evaluate 
this effect under more realistic conditions. We numer- 
ically track the evolutionary history of the universe as 
it shifts from the radiation dominated era to the matter 
dominated era and finally to the cosmological constant 
dominated era. Since the universe is considered to have 
high conductivity after it becomes radiation dominated, 
we make use of the magnetohydrodynamic (MHD) ap- 
proximation [l^l ■ We take into account the Alfven veloc- 
ity, which determines whether or not resonance amplifica- 
tion occurs in the MHD framework fl^ . In addition, nu- 
merical calculations allow us to follow the time evolution 
of the interaction process. We clearly show the behavior 
of the magnetic fields not only in both the long and short 
wavelength limits, but also at the point where each mode 
crosses the horizon, where the impact of the gravitational 
waves is the largest. Because of this property that the 
amplification is most efficient when the mode crosses the 
horizon, we only investigate the effect on magnetic fields 
produced by inflation which is the only mechanism that 
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is able to generate such fields outside the horizon. 

The outline of this paper is as follows. In Sec. we 
first introduce the basics of the covariant formalism which 
is convenient in describing the interaction between grav- 
itational waves and magnetic fields. We then describe 
equations for the Hubble expansion, magnetic fields and 
gravitational waves, which form the base of our inves- 
tigation. In Sec. IIIH we first arrange the equations to 
be applicable for our numerical calculation, and then per- 
form the numerical calculation to evaluate the magnitude 
of the induced magnetic fields, and show and discuss the 
results. We end with our conclusions in Sec. IIVI 



and its evolution is simply given by the Friedmann equa- 
tion, 

H = HoV^ra-^ + ^ma-^ + I^a, (2) 

where we set a{to) = 1 and the subscript "0" denotes 
the present time. The density parameters fir, ^m, 
denote radiation, matter, and the cosmological constant, 
respectively We use flrh'^ = 4.15 x 10"^ and the WMAP 
5 yr results h = 0.705, flm = 0.274, and = 0.726 [H. 

C. Evolution equation for magnetic fields 



II. BASIC EQUATIONS 
A. The covariant formalism 

The interaction between magnetic fields and gravita- 
tional waves has been well studied using the covariant 
approach. We define the 1-1-3 decomposition of space- 
time by choosing the time direction to be along the ve- 
locity vector of the fundamental observer Ua, which sat- 
isfies MqU" = —1. Then the observer's rest frame is de- 
scribed by the projection tensor hat — gab + UaUb- Pro- 
jecting the covariant derivative Va onto the direction of 
Ua and hab, we define the time derivative and the spa- 
tial derivative as Xa-.-b = u'^'^cXa - b and DcXa-.-b = 
hc'^'ha'^ ■ ■ ■ hb^V dXe...f, respectively. The 3-dimensional 
curl is defined as cmXXa = SabcD^X^ , where Sabc is the 
antisymmetric permutation tensor. 

In the covariant formalism, the basic kinematic quanti- 
ties are defined by decomposing the covariant derivative 

of Ua, 



UaUb, 



(1) 



where aab = £'(bWa), ^ab = D[f,Ua], 6 = = D'^Ua 

and iia = u^VbUa are respectively the shear tensor, the 
vorticity tensor, the volume expansion scalar, and the 4- 
acceleration vector. Since we are interested in cosmolog- 
ical magnetic fields under the influence of gravitational 
waves in the expanding universe, we only consider the 
shear aab and the expansion 0, which correspond to the 
tensor perturbation and the Hubble expansion rate, re- 
spectively. Therefore, we neglect the vorticity w"^ and 
acceleration fields u° throughout this paper. 



B. Background equation 

For the evolution equation of the background space- 
time, we assume a spatially-flat Friedmann Robertson- 
Walker universe, for which the metric is given by ds^ — 
—dt^ -f a^{t)5ijdx^dx^ . In this metric space, the volume 
expansion scalar 9 is related to the scale factor a{t) as 
= 3a/a = 3i?, where H is the Hubble expansion rate 



The evolution equation for the magnetic field Ba is 
derived by combining Maxwell's equations jlJl , 



E. 



(a) 



aabE'' - -ehabE'' + curls, - J/,) , (3) 



Bia) = 'JabB" - -QhabB" - curl^ , 
D'^Ba - 0, 



(4) 

(5) 
(6) 



where = —Jau'^ is the charge density and J/^a) = habJ^ 
is the 3-dimensional current. If we neglect the effect of 
the matter components of the universe, pc and J(^a) can 
be set to be zero. This enables us to obtain the evolution 
equation for the magnetic fields by only combining the 
Maxwell equations [g, resulting in a description of 
the propagation of magnetic field waves in vacuum. 

However, the real universe is filled with charged par- 
ticles, and therefore we need to take into account the 
electric current, which is given by Ohm's law. 



J{a) = criEa+eabyB"), 



(7) 



where a denotes the electrical conductivity and Va is the 
plasma 3-velocity. After a large number of charged parti- 
cles are produced at the time of reheating following infla- 
tion, the conductivity of the universe is considered to be 
high enough to assume ct — > oo, which continues until to- 
day even after recombination of the universe Q- In this 
limit, the MHD approximation is appropriate. In this pa- 
per, we also assume instantaneous reheating after infla- 
tion, allowing us to make use of the MHD approximation 
continuously throughout the numerical calculation. 

The evolution equation for magnetic flelds in the MHD 
approximation is derived by using Eq.([71) to replace Ea 
by Ja in the Maxwell equations [121]. In the usual MHD 
limit, charge neutrality is assumed so that pe ~ 0. Ne- 
glecting the second time derivative of Ba and taking the 
limit cr —> CXI, we obtain 



B 



,a,.-^eBa-Sab.e^^^B^D^.d = .abBK (8) 

The plasma 3-velocity in the third term satisfies the Eulcr 
equation, 



V{a) + VbD^Va + ^6Wa + CTabv'' 



= EabcJ^B', (9) 
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where p is the total fluid energy density. Since the dis- 
placement current i?" is negligible compared to the elec- 
tric current J a in the MHD limit, Eq.Q becomes 

J(,) = curls,. (10) 

Here, we have neglected the expansion term in Eq.Q, 
which is considered to be the same order as and also 
the shear term, CTab, which is much smaller than the ex- 
pansion Q. Linearizing and combining Eqs. ([8|), ^ and 
([T0| yields an evolution equation for the magnetic fields 
expressed in terms of only the magnetic field variable B, 
Eqs. (HH), (Ull) and (HH). The de tailed derivation of this 
expression is given in Sec. IIII Al 



D. Evolution equation for gravitational waves 

In most previous work, the evolution of the gravita- 
tional waves is investigated using the evolution equation 
of (Tab H, l9l-[TT|. However, here we choose to use the 
variable hab instead of dab, since we have thoroughly in- 
vestigated the evolution of the gravitational waves using 
hij in our previous work [l5| . This variable is defined 
as the tensor perturbation in a spatially-flat Friedmann 
Robertson- Walker universe, ds^ — —dt^ + a'^{t){5ij + 
hij)dx''dx^ , and its evolution equation is given by the 
Einstein equation, 

h,, + 3Hh,j - -^y^h,, - 0. (11) 

When considering only the transverse traceless compo- 
nents of CTab, which satisfy D°'Uab — 0, hab can be related 
to CTab as (Ta^ = ha^ /2 [1^. Siucc (Tab is a flrst order 
variable, we therefore also treat hab as first order. 

Note that we have neglected the anisotropic stress term 
in Eq. ([TT|) . which is induced by magnetic fields in general 
p7| . This effectively means we assume that the energy of 
the magnetic fields is negligibly small compared to that of 
the gravitational waves. This assumption is valid in the 
initial phase of our calculation, since we consider a sit- 
uation in which the infiation-produced magnetic fields, 
which are generally small, are enhanced by the larger 
gravitational waves. However, after the power of the en- 
hanced magnetic fields becomes comparable to the grav- 
itational waves, the one-sided energy transfer from the 
gravitational waves to the magnetic fields may stop and 
the magnetic fields start to affect the evolution of the 
gravitational waves. Therefore, strictly speaking we need 
to take into account the effect of the anisotropic stress 
term if the energy of the magnetic fields grows to exceed 
that of the gravitational waves. However, since our pri- 
mary aim in the present study is to ascertain whether 
the magnetic fields can be amplified to such a significant 
level at all, (i.e. to the degree suggested by Ref. [6]), we 
are justified in ignoring the anisotropic stress term for 
the time being. 



III. NUMERICAL CALCULATION 

A. Linearization and decomposition of the 
equations 

Here, we rearrange the set of equations, Eqs. dS]), © 
and (flUl) . to a form suitable for numerical analysis, re- 
ferring to the procedure of Ref. [12^. For the numerical 
calculation, we use the linearized and decomposed equa- 
tions. First, we linearize the equation by splitting the 
magnetic fields into a homogeneous background field B, 

which scales as B = — 29i3/3, and the perturbed vari- 
ables. In addition, we split the perturbed magnetic fields 
into parallel Ba and perpendicular B components. 

Ba = B [(1 + B)ea + Ba] , (12) 

where B = {e^Ba - B)/B and Ba = NabB''/B. The 
spacelike unit vector e°, which satisfies Cac" = 1, is par- 
allel to the background magnetic field B, and the projec- 
tion tensor Nab = hab — SaS-b is the metric of the hyper- 
surface orthogonal to e°. In this paper, B is treated as 
first order and B and Ba are second order. Note that the 
perturbed variables of the magnetic fields, B and Ba, are 
normalized by the magnitude of the background magnetic 
field B. This means the linearized equations are valid as 
long as these values are much smaller than 1. We there- 
fore have to be careful that the amplified fields do not 
exceed the background field, B,Ba > 1, in which case 
the linearization method is no longer applicable. Here, 
however, for similar reasons to those stated in Sec. IIID[ 
we apply linearization as a first step since we are primar- 
ily interested in estimating whether the amplification is 
significant.^ 

Similarly, we decompose the shear perturbation into 
three parts as 

aab = ^ab + 2S(aefc) + ^eaCf, - ^Nab^ S, (13) 

where E„fc = {N^a'Ntf - NabN^'' f^Wd, = Nabcr^'^e, 
and E = CTahC^e''. Furthermore we split the spatial 
derivative operator into D\\ = e°'Da = d/dz and D"^ = 
N''^V)b. 

Combining Eqs. ([9|) and ([T0|) and neglecting second 
order terms of Va and Uab, we obtain 

V(a)+\QVa = -{B''DbBa-B''DaBb). (14) 



For the same reasons, we point out that the linearization method 
in Ref. Q is also not applicable if significant amplification truly 
occurs. Furthermore, the effect of the anisotropic stress term 
should be taken into account as discussed in Sec. Ill Dl 
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Taking the time derivative of Eq. ([8]), 

Ba + ^OBa + ^QBa + B^D^Va - BaD%b 

-\QBbD\,a + l^BaD^b + BbD\a - BaD\b 
= CTahS'' + aabB\ (15) 

where we neglect combinations of the spatial curvature TZ 
and Va since they are second order. We then substitute 
Eq. (fT4l) into Eq. ((TSj) and conduct the linearization by 
using Eqs. ([T2l) and (fT3|) and neglecting products of B, 
Ba with first order variables, S, '^ab, Ti- and Va- This 
yields 

(S - ^QB)ea + {Ba - ^eSa) + Cl[-D^Bea 

-{Dpa - D^aD\\B)] + ^Q{ebD''Va - CaD'^Vb) 

= (s-^ei])e, + (i]a-^ei]a), (16) 

where we define the Alfven velocity as C\ = B^ / p. Here, 

we have used B = -295/3 and B = -295/3+4625/9. 
By linearizing Eq. ([5]) , we obtain the equation to rewrite 
the term which includes Va in terms of the magnetic field 
variables and the gravitational wave variables, 

ebD^Va - CaD'^Vb = Bca + Ba - Se, - S,. (17) 

Substituting this and the relation D^^B — —D±aB°', 
which comes from the constraint of Eq. we obtain 
the equations for the parallel component to the direction 
of the background field e°, 

B+^eB-ClD^B^t + ^QE, (18) 

and for the perpendicular component, 

Ba + l&Ba - CliDpa + D^aD^bB') 

= ta + ^e^a- (19) 

Next, we decompose the equations by introducing the 
scalar harmonics Q(k±) which satisfy 

DlQik^)=-^Qik^), (20) 

where k± is the physical wavenumber for the fields per- 
pendicular to e". Taking the divergence and rotation of 
Q(k±)j we define the vector harmonics of even and odd 
parity, 

'^1%.), (21) 



where = eabcs''- These harmonics satisfy -D^Q"^^^ = 
-{kJa)^Qf,^^ and DlQf^^^ = -(kJa^Qf^^y Ve 
then decompose the magnetic field variables in terms of 
the scalar and vector harmonics, 

k± 

- E i^l^M.) + ■ (23) 

k± 

Similarly, the gravitational wave variables are expressed 
as 

k± 

- E(sr^.)^?(^.)+^r^.)^(^.))- (25) 

k± 

We further decompose the variables by introducing the 
scalar harmonics for the parallel component Q{k^^), which 

satisfy Df^Q(^k^^) = ^^^^ ^{k±) 

rewritten as Sf^^^ = Efc,, ^ffc^,fc,|)Q(fc|i)> similarly for 
^ik^V ^Ik^y Hereafter, we drop 

subscripts (fc^,fc||). By applying the above decomposi- 
tions, Eqs. (dHI and result in 

^■^ + ^9S^ + Clir6'5==S^ + ?9I]'^, (26) 
3 3 

similarly for B^^^y and 

B^ + -QB^ + C\4b^ ^t^ + (27) 
3 3 

This is the final expression of the equation for the mag- 
netic fields, which we use in our numerical calculation. 

For the calculation of the evolution of the gravitational 
waves, we follow the procedure of our previous work [l5| . 
We use the Fourier transformed equation of Eq. ((TT)), 

h^ + 3Hh^ + —h^^O, (28) 

where h,j{x) = Ea=+,x / ^''^'^(2'r)"^^^4('^)'*k('^)e''"''- 
The polarization tensors e^'^ satisfy symmetric and 
transverse-traceless condition and are normalized as 
E,,4(£^)* = 25^^'. 

B. Method 

Using Eqs. ©, ^ and we numerically follow 

the evolution of both the magnetic fields and the gravi- 
tational waves under the influence of cosmic expansion. 
Although the magnetic field variables are decomposed 
into several parts, here for simplicity we first follow the 
evolution of the variable B^ . Strictly speaking all field 
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components should be considered. However, since their 
amplitudes are thought to be approximately equal, one 
variable is sufficient to demonstrate the main result of 
this paper - that the final amplification is extremely 
small. The calculation begins at the early stage of the 
radiation dominated universe, from which the MHD ap- 
proximation is valid until the present epoch. The initial 
condition of the induced magnetic field is set to be 
zero for all scales. Note that the variable is the value 
divided by background magnetic field B. This allows 
us to follow the evolution of the magnetic field without 
assuming the strength of the background field. It may 
not be reasonable to assume a homogeneous inflation- 
produced magnetic field as such fields are predicted to 
be a stochastic. However, as mentioned in Sec. IH D| our 
primary aim is to ascertain whether the magnetic fields 
can be amplified or not as suggested in Ref. @. Thus 
following their assumptions, we evaluate this effect un- 
der the simple situation of a homogeneous background 
field. In addition, although the considered situation is a 
little different, as the current term is neglected, the case 
of a stochastic magnetic field is investigated in detail an- 
alytically in Ref. These authors also conclude that 
significant amplification does not occur. 

The initial condition of the gravitational waves is de- 
rived from our numerical calculation presented in Ref. 
[isj . which gives an almost scale-invariant power spec- 
trum. We define the power spectrum of the gravitational 
waves as 



2^ 



E 



(29) 



One may write the power spectrum of the gravitational 
waves generated outside the horizon during infiation as 



'PGW,pTim{k) 



16 / H 



TT \mpi 



(30) 



In our calculation, we assume chaotic inflation, meaning 
that CMB normalization gives Hi^f/nipi ~ 10~^-^. Here, 
the amplitude of the decomposed component is as- 
sumed to contain the entire energy given in Eq. pO|) . 
which means we use = assuming the two polar- 
ization states are equal. We also assume that it prop- 
agates in the direction of the background field, setting 



fcgrav = ^11 in Eqs. (I27p and (j28p . This assumption is in 
fact not strictly true for stochastic inflationary gravita- 
tional waves, meaning that the true amplitude may be a 
factor of a few smaller. However, this does not affect the 
result significantly. 



C. Numerical results and discussion 

We first show the time evolution of the square root of 
the power spectrum in Fig. [1] which is defined as 



Vb 



2^ 



o 




C^=0.5 

gravitational wave 

J I I I L. 





(31) 



FIG. 1: The evolution of B\ ~ yVs shown with the scale 
factor a as the horizontal axis. Here, we show the mode with 
wavenumber fc/ao-ffo ~ 10^. The solid curve shows the case 
where Ca = 0.0, the dashed curve shows Ca = 0.5, and the 
dotted curve shows Ca = 1.0. Additionally, the dot-dashed 
line shows the evolution of the gravitational wave mode hx = 
vPaw- 



This can be interpreted as the power of the induced 
magnetic field perturbation with characteristic scale A — 
27r/fc, Bx sa \^Vb- In this figure, we show the mode with 
wavenumber k — 10^/ao-ffo for three cases of different 
Alfven speeds, Ca = 0.0,0.5, 1.0. Since the energy den- 
sity of the magnetic fields is much smaller than the total 
energy density of the universe, the real universe is con- 
sidered to be close to the case Ca = 0.0. We also plot the 
evolution of the gravitational wave mode h\ w VPgw- 
As seen in the figure, the largest amplification occurs 
at around the point at which a mode enters the hori- 
zon, when the magnetic fields are amplified to the same 
power as the gravitational waves. After the mode en- 
ters the horizon, its behavior depends on the value of the 
Alfven speed. For Ca = 0.0, the induced magnetic field 
maintains a constant amplitude, while for Ca — 0.5 and 
1.0, the induced field undergoes oscillatory decay, with a 
decay rate that depends on the Alfven speed. 

According to Eq. (P7|) . the behavior described above 
can be interpreted as follows. When the mode is outside 
the horizon k ^ aH, the evolution of the magnetic fields 
depends on the amplitude of the gravitational waves as 
-Bk oc Sk. At the same time, the gravitational waves 
remain at constant amplitude when they are outside the 
horizon, Ek ~ /ik ^ 0. This means the magnetic fields 
are barely induced at all by those gravitational waves 
outside the horizon. However, when the mode becomes 
comparable to the horizon scale (fc ^ aH), where the 
gravitational waves begin to oscillate, Ek becomes larger 
and induces magnetic field perturbations. This is the 
reason that the magnetic field perturbation is seen to 
increase toward the time of horizon crossing in Fig. [T] 
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log(k/aoHQ) 

FIG. 2: The square root of the power spectrum of the induced 
magnetic fields B\ = \/Pb at the present time. The horizon- 
tal axis is the wavenumber normalized at the present horizon 
scale (~ 1/aoHo)- The upper panel shows the small Alfven 
speed case, Ca ~ 0.0,0.01,0.1,0.5, and the bottom panel 
shows the large Alfven speed case, Ca = 0.5, 0.9, 0.99, 1.0 



i.e. the case of Ca = 1-0, resonant oscillation occurs and 
pushes up the amplitude of the magnetic field perturba- 
tion. On the other hand, cosmic expansion reduces the 
energy of the oscillating magnetic waves at a rate pro- 
portional to a~^/^ and the gravitational waves at a rate 
proportional to a^^. Since the gravitational waves decay 
faster than the magnetic field perturbations, the amplifi- 
cation effect by resonance oscillation becomes weaker as 
the universe evolves. Therefore, as seen in Fig. [U the 
magnetic field perturbations decay after they enter the 
horizon even when the resonant condition is satisfied. 

Figure [5] shows the square root of the power spectrum 
at the present time for several different values of the 
Alfven speed. We see that the induced magnetic fields 
are the largest in the case of Ca — 0.0. In this case, 
as mentioned above, there is no resonant amplification 
and no decay of the amplitude due to the cosmic expan- 
sion after the horizon crossing, so the amplitude of the 
induced magnetic field only depends on the that of the 
gravitational waves at the horizon crossing. Since the 
primordial power spectrum of the gravitational waves is 
scale-invariant, the power spectrum of the induced mag- 
netic fields also ends up scale-invariant. On the other 
hand, in the other cases where the Alfven speed is not 
zero, the magnetic fields suffer damping due to the cos- 
mic expansion after they enter the horizon. As seen in 
the upper panel of Fig. [21 the damping arises on modes 
smaller than uH/Ca- In the bottom panel, which shows 
the cases where the Alfven speed is large enough to bring 
on the resonance with the gravitational waves, we see 
the amplification pushes up the amplitude of the smaller 
modes. The change in the frequency dependence seen 
around k/aoHo ^ 10^'^ is considered to be due to the 
change of the Hubble expansion rate at matter-radiation 
equality, which causes the change of the frequency depen- 
dence of the amplitude of the source term (gravitational 
waves). In any case, the magnitude of the induced mag- 
netic field is less than 10"^ times the magnitude of the 
original background field at all scales, which leads us to 
conclude that the impact of the gravitational waves is 
negligibly small. 



The evolution inside the horizon k » aH depends on 
the value of the Alfven speed. Setting the Alfven velocity 
to zero is equivalent to saying there is no Alfven mode. 
Thus, the resonant amplification and the decay of oscil- 
lating waves due to cosmic expansion do not occur. This 
is the reason why the evolution of the magnetic field fluc- 
tuations in the case Ca = 0.0 is fiat after entering the 
horizon in Fig. [TJ In this case, the amplitude of the in- 
duced magnetic field is determined only by the amplitude 
of the gravitational waves at their horizon crossing. In 
contrast, if the Alfven mode exists, Ca 0.0, the mag- 
netic field starts to oscillate after entering the horizon. 
In this case, the equation of the magnetic field takes the 
form of a forced oscillator sourced by the gravitational 
wave. In the case where the velocities of the magnetic 
field mode and the gravitational wave mode are the same. 



IV. CONCLUSION 

We have numerically estimated how primordial mag- 
netic fields are affected by inflation-produced gravita- 
tional waves in the framework of the MHD approxima- 
tion. Our calculation has shown that the magnetic fields 
gain most energy from the gravitational waves when each 
mode enters the horizon, although the impact is not large 
enough to call it "amplification". We have found that 
the magnitude of the induced magnetic field is less than 
10^^ times the magnitude of the original background 
field. This conclusion is consistent with the analytical 
estimation in the previous work of Ref. • 

Although the induced magnetic fields are quite small 
in all cases, we have found that the Alfven speed makes 
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a difference to the evolution of the magnetic fields inside 
the horizon. The Alfven speed is the factor which deter- 
mines whether resonance amplification occurs or not, and 
we find that it occurs when the Alfven speed is the same 
as the speed of light. However, this situation would not 
be seen in the real universe, because such a large Alfven 
speed corresponds to a situation in which the magnetic 
fields dominate the energy density of the universe. Fur- 
thermore, even if we assume this unrealistic situation, 
our investigation has found that the amplitude of the 
induced magnetic fields does not increase because both 
the magnetic field waves and the gravitational waves lose 
their energy via the cosmic expansion after they enter the 
horizon. On the other hand, if the Alfven speed is neg- 
ligibly small, the induced magnetic field does not decay, 
which results in a larger amplification rate compared to 
the case of finite Alfven speed. 

One thing we must note is that our calculation does not 
take into account the epoch before reheating, where the 
conductivity of the universe is low, and the epoch dur- 
ing reheating, where the Hubble expansion rate behaves 
as in a matter dominated universe. When the conduc- 
tivity is low, we cannot apply the MHD equations we 
have used in this paper. However, we do not think the 
magnetic fields are significantly affected by the gravita- 
tional waves during the low conductivity epoch, because 
most of the modes we observe today are outside the hori- 
zon during the early epoch, and the gravitational waves 
maintain a constant amplitude. Although the equation 
for magnetic fields in a vacuum is different from that 



of the MHD limit, the source term consists of the same 
components, namely the time derivatives of the gravita- 
tional waves, cr^ = /i^/2 and a\ (see Eq. (3) of Ref. 
Thus, the gravitational waves barely induce the magnetic 
fields outside the horizon, even if the conductivity is low 
in the early epoch. Also, the change in the Hubble ex- 
pansion rate does not affect the behavior of the constant 
gravitational waves when they are outside the horizon. 
Therefore, the early epoch before and during reheating 
is considered to not change our results. 
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